Dynamics of a strongly interacting Fermi gas: the radial quadrupole mode 
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We report on measurements of an elementary surface mode in an ultracold, strongly interacting 
Fermi gas of 6 Li atoms. The radial quadrupole mode allows us to probe hydrodynamic behavior in 
the BEC-BCS crossover without being influenced by changes in the equation of state. We examine 
frequency and damping of this mode, along with its expansion dynamics. In the unitarity limit 
and on the BEC side of the resonance, the observed frequencies agree with standard hydrodynamic 
theory. However, on the BCS side of the crossover, a striking down shift of the oscillation frequency is 
observed in the hydrodynamic regime as a precursor to an abrupt transition to collisionless behavior; 
this indicates coupling of the oscillation to fermionic pairs. 

PACS numbers: 34.50.-s, 05.30.Fk, 39.25.+k, 32.80.Pj 



I. INTRODUCTION 

The advent of ultracold, strongly interacting Fermi 
gases p|, El , molecular Bose- Einstein condensates [3, 0, [E| , 
and fermionic condensates [H, 0] has opened up unique 
possibilities to study the fundamental physics of inter- 
acting fcrmions. The availability of controllable model 
systems with tunable interactions provides unprece- 
dented experimental access to the many-body physics of 
fermionic quantum systems, which is of great fundamen- 
tal importance for various branches of physics [8(. 

A fundamental problem, which has been discussed in 
the theoretical literature for decades 0, EH EH E3> * s 
the crossover from Bose-Einstein condensation (BEC) to 
a macroscopic quantum state in the Bardeen-Cooper- 
Schrieffer (BCS) regime. In this crossover, the nature of 
pairing changes from the formation of bosonic molecules 
by fermionic atoms to pairing supported by many-body 
effects. With novel model systems now available in ul- 
tracold Fermi gases, the BEC-BCS crossover has recently 
stimulated a great deal of interest in both theory and ex- 
periment @. 

Collective excitation modes in trapped ultracold Fermi 
gases provide powerful tools to investigate the macro- 
scopic properties of a system in the BEC-BCS crossover 
pl|. For experiments of this class, ultracold 6 Li gases 
have excellent properties. This is because of their stabil- 
ity in the molecular regime 0, EH EH and precise mag- 
netic tunability of interactions based on a broad Fesh- 
bach resonance 0, E3] • Early experiments on collective 
modes in the BEC-BCS crossover provided evidence for 
superfluidity PH and showed a strikin g tr ansition from 
hydrodynamic to collisionless behavior [l9( . More recent 
experiments yielded a precision test of the equation of 
state [2C| . The previous experiments have focussed on 
collective modes with compression character, where both 
the hydrodynamic properties and the equation of state 
determine the mode frequency PJ, E! EE EH EH El • 

In this Article, we report on measurements of a pure 
surface mode in the BEC-BCS crossover, which provides 
new insight into the dynamics of the system. The "radial 



quadrupole mode" in an elongated trap, the fundamen- 
tals of which are discussed in Sec. HH allows for a test 
of hydrodynamic behavior without being influenced by 
changes in the equation of state. In Sec. IIII1 we present 
our experimental setup and the main procedures. Wc 
introduce a tool to excite collective oscillations with an 
acousto-optic scanning system. The results of our mea- 
surements, presented in Sec. IIV1 provide us with new 
insight on the abrupt transition from hydro dynamic to 
collisionless behavior, first observed in [19j. The present 
work provides strong evidence that quasi-static hydro- 
dynamic theory [24[ does not apply to collective modes 
of a strongly interacting fermionic superfluid, when the 
oscillation frequencies approach the pairing gap [25j . 



II. RADIAL QUADRUPOLE MODE 

The confining potential in our experiments is close to 
the limit of an elongated harmonic trap with cylindrical 
symmetry. In this case, we can consider purely radial 
collective oscillations, neglecting the axial motion. The 
frequencies of the radial modes can be expressed in units 
of the radial trap frequency uj t . We note that our ex- 
periments are performed in a three-dimensional regime, 
where the energy Hu r is typically a factor of 30 below the 
chemical potential and finite-size effects can be neglected. 

In this situation, there are two elementary collective 
modes of the system, the radial compression mode and 
the radial quadrupole mode pH EH- Wc focus on the 
quadrupole mode, which is illustrated in Fig. [T] This 
mode corresponds to an oscillating radial deformation, 
which can be interpreted as a standing surface wave. 
The mode was first demonstrated in atomic BEC exper- 
iments (27j and applied to investigate rotating systems 
[IH, but so far it has not been studied in strongly inter- 
acting Fermi gases. 
Being a pure surface mode, the frequency w q of the ra- 
dial quadrupole mode does not depend on the compress- 
ibility of the system. The frequency uj q does not depend 
on the equation of state but on the collisional proper- 
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FIG. 1: Illustration of the radial quadrupole mode as an el- 
ementary collective excitation of an elongated, trapped atom 
cloud. 

ties. In the hydrodynamic regime, whether the gas is a 
superfluid or a classical gas with a collision rate strongly 
exceeding the radial trap frequency, the frequency of this 
mode is given by [26| 

LU q = V2iU r . (1) 

In contrast, for a collisionless gas, where the atoms freely 
oscillate in the trap, the frequency is 

uj q = 2 u> T . (2) 

Because w q is insensitive to the compressibility of the 
gas and the difference between the collisonless and the 
hydrodynamic frequencies is large, the radial quadrupole 
mode can serve as an excellent tool to probe pure hydro- 
dynamics. Particularly interesting is the transition from 
hydrodynamic to collisionless behavior at lowest temper- 
atures. Such a change occurs in a strongly interacting 
Fermi gas on the BCS side of the resonance [H, HH, HH • 
Near this transition, measurements on the compression 
mode indicated frequency down shifts, which raised ques- 
tions concerning the validity of standard hydrodynamic 
theory in this interaction regime [29|, [3(j. Previous ex- 
periments could not unambiguously identify the origin 
of frequency shifts near the hydrodynamic-to-collisionless 
transition, which is a particular motivation for probing 
the crossover gas with the radial quadrupole mode. 

III. EXPERIMENTAL PROCEDURE 

The apparatus and the basic preparation methods for 
experiments with a strongly interacting Fermi gas of 6 Li 
atoms have been described in our previous work 0, ITjjl 
I2H l3ll |. As a starting point, we produce a molecular 
BEC of 6 Li 2 @, [p. By chan ging an external magnetic 



field, we can control the inter-particle interactions in the 
vicinity of a Feshbach resonance, which is centered at 
834G 0,113 • Thc interactions are characterized by the 
atomic s-wave scattering length a. 

We start our experiments with an ensemble of about 
N = 4 x 10 5 atoms in an almost pure BEC at a mag- 
netic field of 764 G. In order to change the properties 
of the system adiabatically, we slowly ramp to the final 
magnetic field, where the measurements are performed 
[3l| . The temperature of the gas is typically below 0.1 
Tp, unless stated otherwise. 

In order to observe the collective oscillations, we take 
absorption images of the cloud in the x-y-plane after re- 
lease from the trap. We illuminate the atoms with a 
probe beam along the z-direction of the cigar-shaped 
cloud. The probe light causes a resonant excitation of 
the D2-line, at a wavelength of 671nm. We use dichroic 
mirrors for combining and separating the probe and the 
dipole trapping beam. The frequency of the probe beam 
can be tuned over a range of more than 1GHz, which en- 
ables resonant imaging over the whole range of magnetic 
fields that we create in our experiments. 

The gas is confined in a nearly harmonic trapping po- 
tential, which has an axially symmetric, cigar-shaped 
trap geometry. Optical confinement in the radial direc- 
tion is created by a focused 1030-nm near-infrared laser 
beam with a waist of ~ 58 i-j.ni. The potential in the axial 
direction consists of a combination of optical and mag- 
netic confinement [H; the magnetic confinement is dom- 
inant under the conditions of the present experiments. 
We set the laser power to 270 mW, which results in a ra- 
dial trap frequency of u> t ~ 27r x 370 Hz and an axial trap 
frequency of ui z ~2nx 22 Hz at a magnetic field of 764 G. 
The trap frequencies correspond to a Fermi energy of a 
noninteracting cloud E F = h^u^N) 1 ^ = fc B X 740nK. 

In order to excite collective oscillations, we suddenly 
change the optical trapping potential. The position and 
shape of our trapping potential in the x-y-plane can be 
manipulated through the use of a two-dimensional scan- 
ning system. One feature of the system is that we can 
rapidly displace the trap laterally. Fast modulation of 
the beam position enables us to create time-averaged po- 
tentials [33. f33(| . 

The scanning system is constructed by use of two 
acousto-optic modulators (AOMs), which are aligned for 
vertical and horizontal deflections. Fig. [5] illustrates the 
principle of our scanning system for one direction. A col- 
limated beam passes through an AOM and is deflected 
depending on the driving frequency. A lens is placed at a 
distance of one focal length behind the AOM, so that the 
deflection results in a parallel displacement of the beam. 
By changing the driving frequency of the AOM, the lat- 
eral position of the focus is shifted. This system enables 
us to displace the focus of the trapping beam in the hor- 
izontal and the vertical direction by up to four times the 
beam waist in all directions. Furthermore, the deflec- 
tion can be modulated by frequencies of up to ~ 1MHz 
within 3dB bandwidth. In our trap configuration, we use 
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FIG. 2: Schematic illustration of the scanning system. A 
wide collimated beam passes through an AOM. The result- 
ing deflection angle depends on the driving frequency of the 
AOM. The beam passes through a lens at the distance of one 
focal length behind the AOM. The lens focuses the beam for 
atom trapping. A change in deflection angle results in a par- 
allel shift of the beam position in the focal plane. The solid 
and dashed lines show the beam path for different deflection 
angles. The zeroth order beam is not shown. 



modulation frequencies of 100kHz, which greatly exceeds 
the trap frequency. We create elliptic potentials, i.e. po- 
tentials with u; x ^ ujy, by modulating the trap position 
along a specific direction. We use this for the excitation 
of the quadrupole mode. By choosing a suited modu- 
lation function 34], these elliptic potentials are nearly 
harmonic. 

When we excite the quadrupole mode, we first adia- 
batically deform the trapping potential in ~ 100 ms to 
an elliptic shape. This slow deformation ensures that 
the cloud stays in thermal equilibrium even in the near- 
collisionless regime and no excitations occur. We sud- 
denly switch off the deformation leading to an oscilla- 
tion in the x-y-plane of the elliptic cloud in the originally 
round trap. 

The initial deformation corresponds to different trap 
frequencies in horizontal and vertical direction where 
w 0x = (1 + ol)Wt and wo y = (1 — a)to r . The parameter 
a determines the amplitude of the emerging oscillation; 
we choose it for most of our measurements (unless stated 
otherwise) to be a sa 0.05. We increase a by increas- 
ing the modulation for the time averaged potential along 
the y-direction. As the modulation decreases the con- 
finement strength of the dipole trap, we simultaneously 
ramp up the trap power to ensure that the mean trap 
frequency u> T — ^/wox^Oy remains constant. This avoids 
excitation of the compression mode. 

Fig.[3]shows the timing scheme for the excitation of the 
radial quadrupole mode. At t — 0, the collective oscilla- 
tion is excited and the cloud starts oscillating in the trap 
for a variable time itrap ■ Horizontal and vertical widths of 
the cloud, Wx and W y , oscillate in the trap out of phase 
with a relative phase shift of w. As an observable, we 
choose the difference in widths AW = W x — Wy, which 
cancels out small effects of residual compression oscilla- 
tions. For normalization, we introduce the width Wq of 
the cloud in the trap without excitation. 

Experimentally, we determine the collective 
quadrupole oscillations after suddenly switching off 
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FIG. 3: Timing scheme for the excitation of the radial 
quadrupole mode. The ellipticity of the trap is slowly ramped 
up within 100 ms. This results in a change of a in the trap 
frequencies, where a characterizes the ellipticity, and sets the 
initial, normalized deformation AW/Wo = ~2a. Wo is de- 
fined as the width of the cloud in the trap without excitation. 
At t = 0, the elliptic deformation is switched off and the os- 
cillation in the trap begins. (Shown here is an oscillation in 
the hydrodynamic regime.) The oscillation continues until 
the trap is turned off at t — ttrap, which is usually between 
and 10ms. At t = ttrap, the cloud is released from the trap 
and expands for the time £tof , which is typically 2 ms. 



the trap and a subsequent expansion time £tof ■ We 
then take an absorption image of the cloud and deter- 
mine its horizontal and vertical widths W x and W y via 
a two-dimensional Thomas- Fermi profile fit. From these 
measurements after expansion, we can determine the 
in-trap behavior. 

Typical data sets of radial quadrupole oscillations 
are shown in Fig.[4j Fig.rjja) shows an oscillation in 
the hydrodynamic regime; here we observe a weakly 
damped harmonic oscillation centered about a small 
constant offset. Fig.rjjb) shows the typical behavior in 
the collisionless regime. The frequency of the oscillation 
is clearly higher than in the hydrodynamic regime. 
The oscillation shows stronger damping and has an 
exponentially time- varying offset. 

We find that, for both regimes, the dependence of AW 
on ttrap can be well described by the fit function 



AW = A 



COS (Wqttrap 



+ C e c * trap + y , 
which is explained in detail in Appendix [Bj 
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A. Frequency and damping 
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FIG. 4: Typical radial quadrupole oscillations in the hydrody- 
namic (a) and collisionless (b) regimes. The solid lines show 
fits to our data according to Eq.((3]). The dashed lines indicate 
AW = 0. The expansion time £tof is 2 ms. In (a), the os- 
cillation in the unitarity limit (B — 834G) is shown, whereas 
(b) shows the oscillations for B — 1132G (1/fcpa ~ —1.34). 
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Note that the frequency w q and the damping constant 
k are independent of the expansion during £tof and 
characterize the behavior of the trapped oscillating atom 
cloud. In contrast, the amplitude A and the phase shift <p 
depend on the expansion time and provide further infor- 
mation on the dynamics of the gas. The offset function 
C e _, > ttra p with amplitude C and damping constant £ re- 
sults from thermalization effects and is only relevant in 
the collisionless regime (see discussion in Appendix [B| . 
The constant offset tjq results from a slight inhomogeneity 
of the magnetic field, which gives rise to a weak saddle 
potential. This increases (decreases) the cloud size in 
y-direction (x-direction) during expansion. 



IV. EXPERIMENTAL RESULTS 

Here we first discuss our measurements of the fre- 
quency ujq and the damping rate k of the in-trap 
oscillation (Sec. lIV A]) . We then present the data for the 
phase offset <f> and the amplitude A fSec. lIVBj) . Finally, 
we explore the hydrodynamic-to-collisionless transition 
rSec. HVC]) . As commonly used in the field of BEC-BCS 
crossover physics Q , the dimensionless parameter 1 / k?a 
is introduced to characterize the interaction regime. The 
parameter hp = \/2mEp/h is the Fermi wave number 
and m is the mass of an atom. 



FIG. 5: Frequency Lu q (upper plot) and damping rate k 
(lower plot) of the radial quadrupole mode. Both quanti- 
ties are normalized to the radial trap frequency uii and plot- 
ted versus the interaction parameter l/fea. The dashed 
lines indicate the theoretical predictions in the hydrodynamic 
(w q /[j r = v2) and in the collisionless limit (a; q /o;r = 2). The 
shaded area marks the transition from hydrodynamic to col- 
lisionless behavior between 1/fcpa ~ —0.72 (B ~ 930 G) and 
l/fcpa^-0.85 (B^ 960G). 

In Fig.[5j we show the results for the frequency uj q 
and the damping rate k of the radial quadrupole mode 
throughout the BEC-BCS crossover. Both cu q and n are 
normalized to the trap frequency w r , which we determine 
by a sloshing mode measurement |20| . We include small 
corrections resulting from anharmonicity of the trapping 
potential and the residual ellipticity of the trap (see ap- 
pendix [C| . 

The data confirm the expected transition between 
the hydrodynamic and the collisionless regime on the 
BCS side of the resonance (see Sec.|II]). The transi- 
tion is qualitatively different from the hydrodynamic-to- 
collisionless crossover in a classical gas [35f or in a Fermi 
gas without superfluidity [36j]. Instead of a continuous 
and monotonous variation of the frequency between the 
two limits (v^^r and 2w r ), an abrupt change occurs. 
When this transition is approached from the hydrody- 
namic side, a striking frequency downshift shows up as a 
precursor of the transition to higher frequencies. In the 
transition region (shaded area in Fig. [5]), no data points 
are shown because of the large damping and correspond- 
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ingly very large uncertainties for the measured frequency. 

The damping rate shows similar behavior as in our pre- 
vious measurements on the radial compression mode (l9l . 
[23||. Maximum damping occurs near the hydrodynamic- 
to-collisionless transition, whereas minimum damping is 
observed slightly below the resonance. In general, we 
find that damping is roughly two times larger for the 
quadrupole mode than for the compression mode at 
the same temperature |37| . The faster damping of the 
quadrupole mode is plausible in view of the larger fre- 
quency change at the transition. 

We now discuss the behavior in different regions in 
more detail: 

1/fcpa «0: In the unitarity limit, the normalized fre- 
quency agrees well with the theoretically expected 
value of Wq/w r = V2 for a hydrodynamic gas, see 
Eq. (TTJ). To check for consistency with previous 
experiments [2(J, we here also reproduced the fre- 
quency y/10/3 w r of the radial compression mode 
on the 10 -3 accuracy level. The damping is low 
for the Fermi gas in the unitarity limit. In contrast 
to the compression mode, the quadrupole mode fre- 
quency stays constant throughout the crossover, in- 
dicating that it is independent of the equation of 
state. 

l/fcpa>0: In the strongly interacting BEC regime, 
there is an increase in the damping and a slight 
increase in the frequency for increasing 1/fcpa. As 
the gas is more susceptible to heating by inelastic 
processes in the deep molecular regime 13], both 
effects may be due to a thermal component in this 
region. 

1/fcpa w —0.8: The frequency exhibits the pronounced 
"jump" from the hydrodynamic to the collisionlcss 
frequency. This transition is accompanied by a pro- 
nounced maximum of the damping rate. 

1/fcpa < — 0.8: The frequency stays almost constant 
about 5% above the theoretically expected value 
of ujq = 1lo t . Interaction effects in the attractive 
Fermi gas are likely to cause this significant up- 
shift [13, Hil • As we cannot experimentally realize 
a non-interacting Fermi gas above the resonance, 
we could not perform further experimental checks. 

1/fcpa < and 1/fcpa > —0.8 : In this regime, we detect 
a substantial down shift in the quadrupole mode 
frequency. The effect begins to show up already 
slightly above the resonance (1/fcpa = 0) and in- 
creases to a magnitude of almost 20% (w q /o; r » 
1.15 at 1/fcpa = —0.72), before the transition to 
collisionless behavior occurs. Indications of a simi- 
lar down shift have been observed already in com- 
pression mode experiments [HI, [2l|, [23| , but here the 
down shift is considerably larger and not blurred by 
changes in the equation of state. 



A plausible explanation for the curious behavior of the 
collective mode frequency on the BCS side of the reso- 
nance is provided b y c oupling of the oscillation to the 
pairing gap [HI, I25I |29||. If we assume that the abrupt 
transition is caused by pair breaking resulting from reso- 
nant coupling of the oscillation to the gap, then the down 
shift may be interpreted as a coupling effect when the gap 
is not much larger than the oscillation frequency [30(. 
A similar shift may also arise from coupling of hydro- 
dynamics and quasiparticle motion [39]. The observed 
phenomenon still awaits a full theoretical interpretation. 



B. Phase shift and amplitude 

Additional information on the interaction regime is 
provided by the phase shift <f) and the amplitude A of 
the observed oscillation (see Eq.®). This is useful since 
extremely high damping in the transition region makes 
a meaningful determination of frequency and damping 
practically impossible. We find that both amplitude and 
phase shift, however, can be determined with reasonable 
uncertainties even in the transition regime. 

In the following, we present measurements of phase 
shift and amplitude. These are compared to model cal- 
culations, which are described in detail in appendix [A"l 

In Fig. [SI the phase <p an d the relative amplitude are 
plotted versus the interaction parameter 1/fcpa. The rel- 
ative amplitude is given by the amplitude A (definition 
see Eq.([3])) divided by the average width of the cloud 
after expansion. The average width is obtained by aver- 
aging (W x + W y )/2 over one oscillation period using the 
same data set from which we extract A. 

In the transition area around 1/fcpa = —0.8, the phase 
shift <fi shows the step-like change at the transition from 
the hydrodynamic to the collisionless regime. This is 
similar to the jump in frequency in Fig.O In the colli- 
sionless and unitary regimes, the phase agrees with the 
theoretically expected values (solid line and dotted line, 
respectively). 

As a general trend, the relative amplitude is larger in 
the hydrodynamic and smaller in the collisionless regime. 
In the hydrodynamic regime, the relative amplitude de- 
creases for decreasing 1/fcpa, which is explained by the 
change of 7 from 1 to 2/3; 7 is the polytropic index of 
the equation of state (see Appendix [A"|) . At unitarity, the 
relative amplitude agrees well with the numerically calcu- 
lated value for 7 = 2/3 (dotted line). In the collisionless 
limit, the relative amplitude is half of the value at uni- 
tarity, which is also consistent with our calculations in 
App. IA 31 We note that at the transition from the hy- 
drodynamic to the collisionless regime, the value of the 
relative amplitude decreases even below the collisionless 
value. 

In summary, the behavior of the phase shift and the 
amplitude agrees with our model presented in Appendix 
lAl (see also Fig. [H] and Fig. fTU)) . in particular the promi- 
nent change in the phase offset is confirmed. 
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FIG. 6: (a) Phase shift (j> and (b) relative amplitude of the 
quadrupole mode versus interaction parameter l/kva after 
£tof = 2ms expansion. The horizontal lines show calculations 
from our theoretical model: the solid lines in the collisionless 
limit, the dotted lines in the hydrodynamic regime at unitarity 
(7 = 2/3) and the dashed lines in the hydrodynamic regime in 
the BEC limit (7 = 1). These calculated values can be read off 
from Fig. [lO] for the phase and Fig. [9] for the amplitude. The 
shaded area marks the transition between hydrodynamic and 
collisionless behavior between w —0.72 and 1/kFd ~ 

—0.85 (see also Fig.[5j). 



C. Further observations 

The measurements presented in the preceding subsec- 
tions were taken under fixed experimental conditions, 
where only the scattering length a was varied. In this 
subsection we investigate how the transition from hydro- 
dynamic to collisionless behavior depends on the experi- 
mental parameters excitation amplitude, trap depth and 
temperature. 

In a first set of experiments, we explored whether the 
position of the transition depends on the excitation am- 
plitude. We increase or decrease the amplitude by a fac- 
tor of 2. This allows us to compare the oscillations where 
the amplitude is ~ 20 %, ~ 10 % and ~ 5 % of the av- 
eraged width. We do not observe any significant change 
in the position of the transition. 

In general, we find that the transition always occurs 
when the mode frequency is similar to the pairing gap. 
This is supported by the fact that when we vary the 
trap depth the transition occurs at a constant scattering 
length (a w — 5000a , B w 960G) and does not depend 
on 1/fcpa [40( | . A change in laser power of our trapping 



laser influences both Fermi energy Ep and the frequency 
uj q . As we increase the trap power by a factor of 10, we 
also increase the radial trap frequencies by a factor of 
VlO ~ 3.2. This changes the Fermi energy by a factor 
of 2.2 and the pairing gap, which scales like the trap 
frequencies, by roughly a factor of 3 [2S| . These findings 
suggest that the transition is linked to a coupling of the 
collective oscillation to the pairing gap. This is also in 
agreement with earlier results on the radial compression 
mode [Hdl. 
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FIG. 7: (color online) Oscillations of the quadrupole surface 
mode at a magnetic field of 920 G and 1/kpa — —0.66. The 
filled circles correspond to a cold ensemble, whereas the open 
triangles correspond to a heated ensemble. The solid lines are 
fits to the data according to Eq.{3]). 



To explore the temperature dependence of the tran- 
sition between the hydrodynamic and the collisionless 
phase, we use a controlled heating scheme similar to the 
one described in [20j, where we hold the gas in a re- 
compressed trap and let it heat up. We set the mag- 
netic field to 920G (l/k F a = -0.66), i.e. slightly be- 
low the hydrodynamic-to-collisionless transition, where 
the regime is still clearly hydrodynamic. We observe the 
oscillations in a gas at the lowest temperature we can 
achieve in our experiments (filled circles) and in a "hot- 
ter" gas (open triangles) in Fig. [7l The temperature of 
the cold gas is < 0.1 Tf and we believe the temperature of 
the heated gas to be < 0.2 Tp. Figure[7]clearly shows that 
the frequency for the colder ensemble is lower than that 
of the heated one and the amplitude is lower by roughly 
a factor of 2. Using our model in AppendixlAl this indi- 
cates a temperature driven transfer of the ensemble from 
the hydrodynamic to the collisionless regime. 

Thus we find that the radial quadrupole mode is suited 
to detect temperature induced changes of the collisional 
regime of the gas. An exploration of the phase diagram 
of our system depending on temperature is possible, but 
beyond the scope of this article. In our lab, work is cur- 
rently in progress on the radial scissors mode, which turns 
out to be an even better tool for the exploration of tem- 
perature effects. 
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V. CONCLUSIONS 

We have presented measurements on the radial 
quadrupole mode of an ultracold 6 Li Fermi gas in the 
BEC-BCS crossover. As a pure surface excitation, this 
elementary mode probes hydrodynamic behavior with- 
out being affected by changes in the equation of state. 
We have measured the characteristic properties of this 
collective mode in a wide range of interaction strengths. 

Our observations provide new insight into the dy- 
namics of the gas, in particular on the BCS side of 
the crossover, where the character of the oscillations 
abruptly changes from hydrodynamic to collisionless be- 
havior. Our measurements presented in this paper show 
the phenomenon much clearer than in the radial com- 
pression mode [III [U [23[ and provide quantitative data 
on the behavior near the transition. In particular, the 
data show that a substantial down shift of the collective 
mode frequency occurs in the hydrodynamic regime as a 
precursor of the transition. 

The experimental results support the interpretation 
that the coupling of oscillation mode and pairing gap 
[HIIIEIIS] plays a crucial role for the collective excitation 
dynamics on the BCS side of the crossover. We anticipate 
that our new quantitative data on the hydrodynamic-to- 
collisionless transition will stimulate further theoretical 
investigations on this intriguing phenomenon. 
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APPENDIX A: SCALING APPROACH AND 
EXPANSION EFFECTS 

Here we present a theoretical model to describe the 
oscillation of the cloud in the trap as well as its expan- 
sion after release; the model adopts the scaling approach 
applied in [4l|, [4^, |43| ■ The interplay between the dynam- 
ics of the collective mode and the expansion behavior is 
of particular interest as it introduces novel methods to 
investigate the collisional regime. We use a scaling ap- 
proach for both the hydrodynamic and the collisionless 
regime [MIIILII1- In App. lA li the limit of a hydrody- 
namic gas is presented, whereas in App. lA 2l the model in 
the collisionless regime is discussed. Based on these mod- 
els, we show calculated results for the amplitude and the 
phase after expansion in App. lA 3"1 



The scaling approach describes the cloud at the time 
t after excitation [4l|, [42|, [43| . Using the scaling function 
bi(t) for i — x,y, the width Wi(t) for all times t > can 
be written as 



Wi{t) = bi{t)wM, 



(Al) 



where W x (0) = (1 - a)W and W y (0) = (1 + a)W are 
the initial widths at excitation and Wq is the width of 
the cloud without excitation. The initial conditions for 
the scaling function are bi(0) = 1 and 6;(0) = 0. 



1. Dynamic behavior in the hydrodynamic limit 

In the hydrodynamic limit, the equations of hydrody- 
namics lead to the following differential equations for b x 
and by [H 



b x {b x by) 
by (b x b y ) 



> 2 

7 "v w y 



(A2) 



where 7 is the polytropic index of the equation of state 
and b z (t) — 1 for our elongated trap geometry. The pa- 
rameters uox and u>o y are the trap frequencies at the mo- 
ment of excitation (t = 0), when the cloud has no further 
excitation and is in thermal equilibrium. In contrast to 
this, Lu x (t) and uj y (t) are the time dependent trap fre- 
quencies. The timing scheme is illustrated in Fig. [31 The 
following equation summarizes the behavior of the trap 
frequencies u>i(t): 



U>i{t) 



t = 
< t < t 



trap 



(A3) 



t > t 



trap ■ 



This enables us to calculate the scaling functions b x and 
by as solutions of Eq. (|A2|) for the in-trap oscillation. In 
the limit of small amplitudes (a < 1) solutions are 



1 + a(l — cos co q t), 
1 — a{\ — cos Wqt), 



(A4) 



where u> q = ^/2uj t is the radial quadrupole oscillation fre- 
quency. Together with Eq. ([Al~|) . we are able to determine 
the difference in widths of the cloud to be 



AW = -2 aW () cos uj n t. 



(A5) 



2. Dynamic behavior in the collisionless limit 

In the collisionless limit, the following set of two uncou- 
pled equations characterizes bi, where i stands for x,y, 



bf 



(A6) 
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In the limit of small amplitudes (a<l) solutions of the 
in-trap oscillation are 

b x = 1 + "(I - COSLU q t), 
(X 

by = 1- -(1-COSWqf), (A7) 

where uj q = 2w r is the radial quadrupole oscillation fre- 
quency. Together with Eq. (|A1|) , we are able to determine 
the difference in widths of the cloud to be 

AW = -aW (1 + cos u q t) . (A8) 

In contrast to the hydrodynamic limit, the oscillation is 
initially not centered around AW = 0. Furthermore the 
oscillation has an amplitude 1/2 of the amplitude in the 
hydrodynamic gas. 

Besides the finding of analytical solutions, it is enlight- 
en P A 

I 1 x,y 




the dynamics of any point in phase space is a simple cir- 
cular rotation about the origin with frequency uj r . Thus, 
the solid circle in Fig. [5J (a) indicates an equilibrium phase 
space contour for the round trap. Right after applying 
the excitation scheme as described in Sec. lIIII the phase 
space contours in the x- and y- direction are given by the 
dashed and dotted ellipses in Fig. [5] (a). Since the gas 
is fully thermalized at the instant of excitation, the ini- 
tial momentum distribution in x- and y- direction is the 
same. As time progresses, the elliptic contours will rotate 
with frequency oj t (see Fig.[SJ(b)), which corresponds to 
oscillations in the trap. We note that both the spatial 
and the momentum distribution in the x-direction are 
never larger than the ones in the y-direction. Therefore, 
AW oscillates between 2aWo and zero and the aspect 
ratio of the cloud never inverts. This is to be compared 
to the hydrodynamic case where AW oscillates between 
±2aW . 

Residual thermalization effects in a near collisionless 
gas will damp out the initial oscillation amplitude of aWo 
and one will eventually end up again with a circular phase 
space contour (see Fig. [8] (c)). This is studied in detail 
in AppendixlBl 




oo r t = n/2 



CO r t — > 00 

(co q = 2co r ) 



FIG. 8: Phase space dynamics for the quadrupole mode in 
the collisionless regime. Shown are phase space contours of 
an ensemble of particles which is held in a round trap (i.e. 
uj x — u) y = uj r ). In (a) and (b) the situation during the oscil- 
lation in the trap is shown for two different times t. The solid 
line indicates the equilibrium phase space contour (without 
excitation), whereas the dotted (dashed) line shows the con- 
tour in the x (y) direction after excitation of the oscillation 
mode, (c) After long times, residual thermalization finally 
damps out the oscillations and leads to a circular phase space 
contour. 

ning to understand the collective oscillations in the colli- 
sionless limit by considering the phase space dynamics of 
the cloud. In Fig. [5J we show the contours of phase space 
distributions in the x- and y- directions. The axes are 
scaled such that for the round trap, i.e. lu x — uj y = cu r , 



3. Amplitude and phase on expansion 

Here we present our calculated results based on the 
models in App. lA ll and App. lA 21 for the hydrodynamic 
and the collisionless limit, respectively. We show the 
relative amplitude that is given by the amplitude A 
(definition see Eq.([3|)) divided by the average width 
of the cloud after expansion (for definition details see 
Sec. lIVBj) . Calculations of this relative amplitude are 
shown in Fig. [5J whereas calculations and measurements 
for the phase offset <\> are shown in Fig.flTJl 

Fig. [5] shows the calculated relative amplitude of a 
surface mode oscillation in the hydrodynamic (dashed 
and dotted curves) and in the collisionless (solid curve) 
regime as function of the reduced expansion time 
WriTOF- The hydrodynamic curves are calculated for 
the BEC limit of 7 = 1 (upper, dashed curve) and in 
the unitarity limit of 7 = 2/3 (lower, dotted curve). The 
amplitude in the collisionless regime is smaller than in 
the hydrodynamic regime. Initially the amplitude of the 
excitation is half as large in the collisionless as in the 
hydrodynamic regime, as already explained in App. lA 2"! 
In expansion the normalized amplitude stays constant in 
the collisionless regime and in the hydrodynamic regime 
for 7 = 1. For 7 = 2/3 in the hydrodynamic regime it 
decreases for longer expansion times. 

In Fig. [TU]we compare experimental data for the phase 
shift cf> with numerical simulations. The data have been 
taken at unitarity where 1/k-pa = (hydrodynamic, 
open circles), and on the BCS-side of the resonance at 
1/kpa = —1.34 (collisionless, closed triangles). The 
dashed line is based on a model for the hydrodynamic 
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FIG. 9: (color online). Calculated relative amplitude of a 
surface mode oscillation versus reduced time of flight w r tTOF 
after release from the trapping potential. The values are cal- 
culated for the hydrodynamic (dashed curve: 7 = 1, dotted 
curve: 7 = 2/3) and collisionless regime (solid curve). The 
vertical dotted line marks the typical expansion time in our 
experiments. 



FIG. 10: (color online) Phase <j> of the collective surface 
mode as detected by fits according to Eq.© versus reduced 
expansion time u T tTo~F at unitarity (open circles) and at 
l/k^a = —1.34 (filled triangles). The lines are numerical sim- 
ulations for the hydrodynamic (dashed line) and collisionless 
regime (solid line). The vertical dotted line marks the typical 
expansion time in our experiments. 



interaction regime and the solid line on a model for the 
collisionless regime. The data agree with the theoreti- 
cal model where no free fit parameters are used. This 
confirms our approach presented above. 



APPENDIX B: THERMALIZATION EFFECTS IN 
A NEAR-COLLISIONLESS GAS 

Here wc describe thcrmalization effects in a near- 
collisionless gas that are not included in the model for 
the collisionless limit in App. lA 21 Despite the word "col- 
lisionless" , collisions play a crucial role for thermalization 
for our experimental parameters. A typical time scale for 
thermalization processes is only a few oscillation cycles 
long. By analyzing the theory, we are able to introduce 
a universal fit function, as given by Eq.©, which de- 
scribes the oscillation both in the hydrodynamic and in 
the near-collisionlcss regime. 

The measured behavior of the nearly collisionless 
quadrupole oscillation (see Fig. Q| has two characteris- 
tics: after excitation the oscillation is centered around 
AW = (W x (0) - Wy(0))/2, then after some time it is 
centered around AW = 0. These two limits are consis- 
tent with thcrmalization of the gas on a relevant time 
scale greater than the period of the oscillation. 

In order to model these effects, wc follow a theory 
based on a classical gas in the transition between the 
hydrodynamic and the collisionless behavior described 
in [38|. An application of this theory for the compres- 



sion mode in the hydrodynamic regime has been used in 
[4lj . Here we will handle thermalization effects of the 
quadrupole mode in the near-collisionless regime. 

Using the classical Boltzmann-Vlasov kinetic equation 
in the relaxation-time approximation and ignoring mean 
field effects one can derive the following coupled differ- 
ential equations [H[ 

&i = wg,^-«?&i (Bl) 

and 

e i = —{e i -6)-2 b ^e i . (B2) 

tr Oi 

The parameter bi is the scaling function described 
earlier in Appendix [Aj 9i is a scaling parameter directly 
related to the temperature and 9 — |£fct9fc. The initial 
condition for Oi is 6*^(0) = 1, as long as the gas is in 
thermal equilibrium at the moment of the excitation. 
The parameter tr is the relaxation time which describes 
the time scale of collisions. In the collisionless limit, 
when tr — > 00, the differential equations (|B1[) simplify 
to the simple form in Eq. (|A6|I . For the hydrodynamic 
limit (tr -> 0), we find Eq. (|A2|) . 

The solutions to these equations depend on the param- 
eter tr as can be seen in Fig.[TTJ Our measured data in 
the collisionless regime are well described by lu^tr ~ 2.3 
(compare to Fig. d]). 
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FIG. 11: Calculated quadrupole oscillations in the near- 
collisionless regime. The lines show the relative difference 
in width AW as a function of the reduced time a> r itra P . The 
oscillation is modeled according to Eq. ljBip and (|B2[) . The 
dark line shows the result of the calculation when lo t tr — 2.3 
and the grey line shows the coscillation in the collisionless 
limit at lj t tr — 1000. 

The universal fit function 

We find that the model calculations from (|B1[) and 
(|B2|) can be well described with the following fit function 

AW = Ae~ Ktt '^ COS (Wqttrap + <t>) 

+ C er^v + y . 

The first term describes the exponentially damped oscil- 
lations. The second term describes the shift of the center 
of the oscillation in the collisionless regime. The third 
term yo is a constant offset which will be discussed later. 

We have used Eq. (|B3|) to fit our experimental mea- 
surements. We find that the free fit parameters £ and k 
are related through £/k « 1.5 for all our measurements 
in the near-collisionless regime. In the hydrodynamic 
regime C = 0, and therefore £ becomes irrelevant. 

The constant offset yo is due to an experimental 
artifact that results from a slight inhomogeneity of 
the magnetic field. At the location of the atoms the 
inhomogeneous magnetic field leads to a weak saddle 
potential which causes a slight anisotropic expansion 
during time of flight. This anisotropy is responsible for 



a slight offset in AW. 

APPENDIX C: CORRECTIONS TO THE 
NORMALIZED FREQUENCY 

The theoretical normalized frequencies Lo q /u) r are cal- 
culated for perfectly harmonic trapping potentials in an 
idealized symmetric trap geometry. There are small 
derivations from this conditions in real experiments. In 
order to compare the experimental data to the ide- 
alized theoretical case, we have to correct our data. 
The measured normalized frequency uj q /u) T of the radial 
quadrupole mode has to be increased because of two 
small corrections. The larger correction is based upon 
a slight anharmonicity of the trapping potential and the 
spatial extension of the cloud in the trap. The smaller 
correction is caused by a small residual ellipticity of the 
trapping potential. 

The potential created by our trapping beam has a 
Gaussian shape. This results in a nearly harmonic po- 
tential in the center of the trap; however, for higher pre- 
cision one must take into account higher order terms of 
the potential. Anharmonicity effects influence both our 
measurements of the sloshing mode frequency, where we 
determine w r , and our measurements of the quadrupole 
mode frequency uj q . As we evaluate the normalized fre- 
quency uj q /Lu r , the anharmonicity effects on sloshing and 
quadrupole mode almost cancel out each other. The 
small remaining correction to the normalized frequ ency 
is included by multiplying with a prefactor 1 + ba [2fJ,[4J|. 
The anharmonicity parameter a relates the energy of the 
oscillation to the total potential depth and is defined by 
a = imw 2 r 2 ms /Vb, where r rms is the root-mean-square 
radius of the trapped cloud and Vq is the potential depth. 
The parameter b depends on the interaction regime. In 
the hydrodynamic regime, it is given by (4+107)/ (2+77), 
whereas in the collisionless regime b is determined by 6/5 
[44| . Here, 7 is the polytropic index of the equation of 
state. In our experiments, typically ba « 0.014 , but bo 
can rise to an upper limit of ba < 0.027. 

In the hydrodynamic regime, there is also a correction 
due to residual ellipticity effects. This correction takes 
into account that we compare our measurements with a 
theory for non-elliptic geometries. The ellipticity e of the 
trap is defined by e = (u> y — u> x )/uj r . In our experiments, 
the ellipticity is small and given by e w 0.07. Therefore, 
we can apply the ellipticitycorrection by multiplication 
of a prefactor 1 + Ae 2 [2(1 0] , where the interaction de- 
pendent factor A is given by (7 + 2)/(47). Altogether, 
Ae 2 is smaller than 0.006 for all data points. 
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